setwd(wkdir)

today <- '20210211'
##
## FISHNET
##
fishnet_lake_chad_ntl_x_d01_shp <- paste0(wkdir,'/local/gis_data/003_boundaries/fishnet/fishnet_lake_chad_ntl_x_d01')
fishnet_only <- read_sf(dirname(fishnet_lake_chad_ntl_x_d01_shp),basename(fishnet_lake_chad_ntl_x_d01_shp))

##
## NTL ZONES
##
ntl.zones <- read_sf(dsn=paste0(wkdir,'/local/gis_data/003_boundaries/fishnet'),layer='fishnet_lake_chad_ntl_adm0_x_d01')
countries_GID_0 <- unique(ntl.zones$GID_0)
# add area km2 #
ntl.zones$sf_area_km2 <- st_area(ntl.zones) / 1000000

out.df <- as.data.frame(st_drop_geometry(subset(ntl.zones,select=c("ID","OBJECTID","GID_0"))))
summary(out.df)
save.dta13(out.df,paste0(wkdir,'/proc_data/CROSSWALK',today,'.dta'))
